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Abstract 

An analytic formula describing the suppression of the equilibrium melting line by quenched point 
pinning disorder is compared to data from ac susceptibility and magnetization measurements in the 
mixed phase of the layered dichalcogenide low-Tc superconductor 2H-NbSe2. This material exhibits 
a sharp peak effect in the critical current jc close to the upper critical field Hc2(T). Arguing that the 
disorder-suppressed melting line in this system is to be identified with the locus of peak positions of 
the critical current as magnetic field and temperature are varied, we demonstrate that this formula 
provides a remarkably accurate fit to the experimental data over three orders of magnitude in 
magnetic field. 
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The translational order of the Abrikosov flux hne lattice is disrupted both by thermal 
fluctuations and by quenched disorder [THl]. Thermal fluctuations drive a flrst-order flux- 
lattice melting transition in pure systems while quenched disorder destabilizes translational 
long-range order even at zero temperature [HI E]. A delicate balance between elastic restor- 
ing forces and the twin disordering effects of thermal fluctuations and random pinning thus 
determines flux line structure and phase behavior in the mixed phase of type-II supercon- 
ductors with quenched disorder [3]. 

Translational correlations cannot be truly long-ranged in a crystal with quenched 
disorder[6]. Qualitatively different states may, however, be identifled through the nature of 
the decay of such correlations P-[3]. In the Bragg glass (BrG), a quasi-lattice phase believed 
to be stable at weak disorder and low temperatures, such correlations decay asymptotically 
as power laws [71 |H]. Upon increasing either the effective disorder strength or the amplitude 
of thermal fluctuations, the Bragg glass transforms into a phase in which such correlations 
decay exponentially. The disordered liquid (DL) is one such stable state obtained at large 
temperatures (T). Another disordered state is obtained in fleld (H) scans out of the Bragg 
glass at low T[3]. Such a disordered state has glassy attributes, including divergent time- 
scales for structural relaxation and relatively short-ranged correlations. It may represent a 
new thermodynamic phase, the vortex glass phase, distinct from the equilibrium disordered 
liquid[21[Sl[in]. 

At large H, the DL phase appears to transform continuously into the glass when T is 
reduced. At low and intermediate values of H, it appears to freeze discontinuously. This 
freezing in the intermediate interaction-dominated regime has traditionally been understood 
as occurring directly into the Bragg glass[lll |12]. Such a view is, however, at odds with 
a large body of data in which anomalies associated with the freezing transition, such as 
a discontinuity in the magnetization, often appear well separated from another structural 
transition which occurs at lower temperatures, particularly for more disordered samples [T3l- 
[T3] . An attractive way to reconcile a large number of experimental observations on both 
high-Tcand low Tc materials is via the proposal that two transition lines always separate 
the BrG phase from the DL phase, with a sliver of disordered glassy phase intervening. 
This possibility is illustrated in Fig. [l |l3f[T5] . with the intermediate glassy phase termed 
as a "multi-domain glass" . The phase diagram of Fig. [T] has been argued elsewhere to be 
a generic phase diagram for the mixed phase with quenched point pinning [T5]. Among its 
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distinctive features is the identification of a smooth connection between high-field and low- 
field glassy phases, a proposed equivalence of the peak effect as seen in temperature scans 
as T — )■ Tc and the "fish-tail" effect seen in field scans at low temperatures, as well as a 
specific prediction for vortex-line structure in the intermediate glassy phase [T314T5] . 

In Fig. [T| the Bragg glass melts into the disordered liquid phase on increasing T via an 
intermediate glassy phase [T3HT5]. This glassy regime broadens both at large H as well as 
very small if, reflecting the increased importance of disorder both at high fields, where the 
multi-domain glass phase is encountered, and at low fields, where a reentrant disordered 
state has been predicted and indeed seen [16]. The intervening sliver of disordered phase 
has been identified with the peak effect regime, the narrow region close to Hc2 in (if, T) 
space in which the critical current jc increases anomalously in a variety of low and high-Tc 
superconductors [T71 [T91422] . Susceptibility measurements on a variety of low-T^ materials 
indicate that the phenomenon of the peak effect itself exhibits a generic two-step character, 
in agreement with the proposed phase diagram. In this picture, the first phase boundary 
encountered when cooling from the disordered liquid state is the remnant of the thermal 
melting line in the pure system, renormalized suitably by disorder. This phase boundary is 
then identified with the loci of peaks in jc in (if, T) space. 

Significantly, this proposal also motivates a new interpretation of the classic problem of 
the peak effect: the peak effect is simply a result of the disordered, high jc vortex glass phase 
intervening between two relatively low jc phases, the Bragg glass phase and the equilibrium 
disordered liquid. In this picture, the "hardening" of the vortex system as a consequence 
of the transition into the multi-domain glass phase (rather than the smooth "softening" 
envisaged in Pippard's original scenario[2I?J or a possibly sharp crossover from collective to 
individual pinning [2U [25]), is the underlying cause of the peak effect [H]. Such an increase in 
critical currents in poly-domain structures has indeed been seen in recent simulations [27]. A 
recent study of this intermediate phase is reported in Ref . |28] , while earlier work examined 
this regime using scanning ac Hall microscopy, proposing the possibility of "coexistence" 
within the peak effect regime as an explanation for the anomalies seen within this regime [29]. 

While translational correlations in the BrG phase are of power-law character, those in the 
DL phase are short-ranged, with typical scales of order a few interline spacings. We have 
argued elsewhere that the most appropriate description of structure in the intermediate 
disordered phase is in terms of a disordered arrangement of ordered domains, as in a "multi- 
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Temperature 

FIG. 1. [Color Online] Generic phase diagram for type-II superconductors with point pinning 
disorder. The Bragg glass is depicted as melting directly into a "multi-domain" glass phase and 
only then into the disordered liquid upon increasing the temperature. The Hp;(T) (where pi stands 
for "plastic flow" ) and Hp lines, associated with the onset and the maximum of the peak effect 
anomaly in the ac susceptibility or critical current, are identified with the transitions between 
the Bragg glass and the multi-domain glass and the multi-domain glass and the disordered liquid 
respectively. The boundary between Bragg Glass and multi-domain glass phase is reentrant [5l 116]. 
although the regime of reentrance appears to be very small in experiments (it is exaggerated for 
clarity). The regime of intermediate vortex glass phase is also exaggerated in the figure vis a vis 
the experiments, where it appears as a very narrow sliver in the phase diagram, concomitant with 
other anomalies signaling a sequence of equilibrium phase transitions [13]. 
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FIG. 2. [Color Online] Experimental data (filled circles) showing the locus of points associated with 
the peak of the peak effect anomaly, as obtained from ac susceptibility measurements. The solid 
line passing through these point corresponds to the theoretical prediction in Eq. [5]for the disorder- 
induced suppression of the phase boundary between the multi-domain glass and the disordered 
liquid. The parameter values used in the fit are C = 10^, m = 10.0, Tc = 7.09 and a = 1.33, 
assuming XY-type critical behaviour. The dashed line represents a theoretical fit which corresponds 
to the pure system behavior, in which we set m = 0. 

domain glass" (MG) phase [T314T5] . Using formulae appropriate to the collective pinning 
regime, a conjecture for typical domain sizes, and experimentally obtained values of jc 
in very weakly disordered samples of 2H-NbSe2 indicate Rd/a ~ 10^, with a the mean- 
intervortex spacing, suggesting that the domain sizes in the MG phase of systems with 
low levels of pinning can be far larger than the correlation length at freezing in the pure 
system[T3l E]. Such a picture rationalizes the association of magnetization discontinuities 
with the transition out of the DL phase. Further, it suggests that examining the instability 
of the DL phase to a phase with solid-like structural correlations may be a good starting 
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FIG. 3. Real part of the ac susceptibility of an NbSe2 crystal measured at three representative 
fields applied parallel to the c-axis. The data is measured using an ac field amplitude of 1 Oe and 
10 Hz frequency. 

point for the calculation of the DL-MG phase boundary. 

This paper compares a generalized version of the predictions of a simple semi-analytic 
theory of the DL-MG phase transition, proposed in Ref. [H], with results from ac sus- 
ceptibility and magnetization measurements on the low-Tc dichalcogenide superconductor 
2H-NbSe2. This material, with Tc — 7. IK, shows a remarkably sharp peak effect signal in 
an interval within about 10% of Hc2(T). Structure in the peak effect regime in 2H-NbSe2 
has been conclusively demonstrated to be domain-like, validating the approach here |30j : for 
related simulations see Ref. p7] . 

We find that a relatively simple analytic formula provides a remarkably accurate parametriza- 
tion of the MG-DL transition line, the line in the phase diagram of Fig. [T| provided, 
as argued extensively elsewhere, the loci of the peak in jc is identified with this transition. 
This fit is shown in Fig. |2| using data from ac susceptibility measurements. Fig. |3] exhibits 
representative data (see below) for the ac susceptibility, while magnetization-based mea- 
surements of the critical current are shown in Fig. |4] at one value of the temperature. As 
can be seen from Fig. [2} the fit to the data is of exceptional quality. 

The approach of Ref.[ll], described briefly below for completeness, draws from a early 
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FIG. 4. Critical current of the 2H-NbSe2 sample as a function of applied field as obtained from 
dc magnetization measurements, illustrating the peak effect in this material, at a temperature of 
T = 5.8K. 



replica theory of correlations in disordered fluids and a replica generalization of the den- 
sity functional theory of freezing as applied originally to the pancake vortex system in 
BSCC0[2Il Such density functional theories have been shown to provide a quantita- 
tively accurate picture of the freezing transition in the pure system, including such details as 
the physics of the anomalous slope of the melting line in these systems and a rationalization 
of asymmetric hysteresis across the fiux-lattice melting transition as induced by the presence 
of free surfaces [STj [38l 139] . 

Replica methods are generalized to disordered fluids through the model problem of 
a system of point particles interacting via a two-body interaction and an explicit one- 
body disorder term Vrf(r)[32]. Here Vrf(r) represents a quenched, random, one-body po- 
tential, drawn from a Gaussian distribution of zero mean and short ranged correlations: 
[Vc;(r)Vrf(r')] = K{\ r — r' |), with [■ ■ ■] denoting an average over the disorder. 

Two point correlation functions, within replica theory, carry two replica indices, denoted 
by the Greek letters a and /3. Thus, the pair correlation functions h°'^{r) and the re- 
lated direct correlation functions C"^(r) characterize the replicated system at the level of 



two-particle correlations. Assuming replica symmetry, as appropriate to the equilibrated 
disordered fluid, we define the replica off-diagonal (a 7^ (3) and replica diagonal (a = (3) 
correlation functions, writing C"/^ = C^^^d^^ + C^^^l-d^iB) and /i"^ = h^S^^ + h^'^^l-d^ii). 
Physically, the function h^^^ describes the disorder-averaged equal-time (equilibrium) cor- 
relation of fiuctuations of the local density, whereas h^"^^ represents the disorder-averaged 
correlation of disorder-induced deviations of the time-averaged local density from its aver- 
age value pe- Direct simulation tests of the predictions of the theoretical framework for these 
correlation functions are reported in Ref. 

The density functional theory of freezing takes such correlations as input to a mean-field 
determination of when the fiuid, with density pe becomes first unstable to the formation 
of a static density inhomogeneity. Applying replica analysis leads to a density functional 
which qualitatively resembles the density functional of a pure system, but with renormalized 
correlations! 



dr 



p(r) In ^ - 6p{r) 



-\jdvj dv'[C^^\\ r - r' I) - C^'\\ r - r' |)][p(r) - p,][p(r') - p,] + . . . . (1) 

Assuming p'^(r) = p(r) for all a, the properties of the density functional are governed by an 
effective direct correlation function given by C^^f{\ r — r' |) = C7''"'^''(| r — r' |) — C^'^\\ r — v' |). 
In this description, disorder enters through (i) the suppression of r — r' |) and (ii) the 

non-trivial character of C'-^-'d r — r' |). 

The inter-replica interaction PK{p,nd) appropriate to vortex lines and pancake vor- 
tices is obtained assuming the principal source of disorder to be atomic scale pinning cen- 
ters [3S]. A model calculation yields /3V^'^\p,nd) = —(3K{p,nd) ~ — F' exp(— p^/^^)5„_o, 
where l3V^'^\p,nd) = l3V°'^{p,nd) with a ^ fi. — IsA is the coherence length in the ah 
plane, and F' ^ lO^^F^ for point pinning of strength drlH^/Sir, with d the interlayer spac- 
ing (~ 15 A) for a layered superconductor (we set = ^ for an isotropic superconductor), 
F = l3d^l/4:TcX^ and (3 = l/ksT. Defect densities of the order of 10^°/cm^ are assumed [32]. 
We will use only the dependence of F' on the field and temperature in our discussion i.e. 
F' ~ F2 ~ 1/T2. 

Further progress requires a calculation of the correlation functions C^^\r) and C^'^\r). 
A replica generalization of the Ornstein-Zernike equations coupling C^^\r) and C*^^-*(r) to 
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h^^^)(r) and h^'^\r) can be derived: these equations must be supplemented, as usual, with 
appropriate closure schemes, such as the hypernetted chain (HNC) or Percus Yevick (PY)[32l 

Eg. 

In mean-field theory, the freezing transition of the pure system occurs when the density 
functional supports periodic solutions with a free energy lower than that of the uniform 
fiuid[33J. We now specialize to the fiux-line lattice, in which case properties of the freezing 
transition are controlled by a two-dimensional C'^-^-^{q±,qz = 0), the Fourier transform of 
C^-^-^(r_L, z). Setting to zero simply corresponds to considering arrays of pancake vortices 
in perfect registry or, equivalently, straight vortex lines representing the crystal [31] . 

The Hansen- Ver let criterion for the freezing of a two-dimensional liquid [351 EZj indicates 
that freezing occurs when the structure factor S{q) = 1/(1 — pC^^^ {q)), evaluated at g.^ i-e. 
S{qm), attains a value of about 5, a value roughly independent of the nature of the inter- 
action potential. A full density functional calculation yields the same quasi-universality, a 
direct consequence of the fact that correlations, not microscopic potentials, are the principal 
determinants of freezing[33]. Since C^'^\qm) > 0, and C^^\qm) is always reduced (although 
weakly) in the presence of disorder, the equilibrium melting fine is always suppressed by 
quenched disorder. 

We use the following ideas: (i) The diagonal direct correlation function C^^\q) is weakly 
affected by disorder and can thus be approximated by its value in the absence of disorder, 
(ii) The off-diagonal direct correlation function C*-^-'(g) varies strongly with disorder and 
with H and, (iii) C^^\q) is well approximated at g = = ^n/a by its value at g = 0. Since 
the Hansen- Verlet condition is satisfied along the melting line, the following holds: 

p.C^" = peiC^'Kqm) - C^'\qm)) ^ 0.8. (2) 

The off-diagonal correlation function C^^-* decays sharply in in real space; in Fourier space, 
therefore, its value at g = g^ is close to its value at g = 0. The vahdity of such an 
approximation has been tested in Ref. plj. Thus, 

C(2)(^) ~ -/3V^^)(r). (3) 

The prefactor scales with temperature as F^, implying that piC^'^\qm) ~ Note that 
C*-^''(gm) increases as B is increased or as T is decreased, as is intuitively reasonable. 

We observe that C^^\qm) increases with a decrease in T; reducing T increases correlations. 
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The variation in C*-^^ is expected to be smooth within the equihbrium disordered hquid. Thus 

C(i)(g„; T - AT, B) = C^^\qm-, T, B) + T) AT, (4) 

where q{B,T) [q > 0) is a smooth function of B and T close to the melting line. To 
first order in C^'^\ {B,T) can be replaced by {Bm,,Tm) and C'^^qm] Bm,Tm) by its value at 
freezing for the pure system: pC^^\qm', Bm,Tm) — 0.8. We may also neglect the B and T 
dependence of q i.e. q{Bm,Tm) — q, with q a constant; at melting this dependence should 
be small provided a ^ A. 

An approximate expression for the suppression of the melting line from its value for the 
pure case {Bm,Tm) = {BmiT),T) now follows: AT^ ~ ^"^P = m '^'^2"^^ , with m a constant, 
related to a temperature derivative of the direct correlation function. Here AT^ is the 
shift in the melting temperature induced by the disorder. This relation predicts a larger 
suppression of the melting line at higher fields and lower temperatures. 

The above result parametrizes the suppression of the melting line by quenched disorder. 
This result can be combined with results from a calculation of the melting line in the pure 
system to obtain a simple analytic formula for the MG-DL phase boundary [H]. At low 
fields, a simple Lindemann parameter-based calculation of this phase boundary appears to 
be reasonably accurate and yields Bm(T) = C(T — Tc)"", where Tc is the critical temperature 
and C is a constant appropriate to the pure systempQ. The exponent a is characteristic of the 
fluctuation regime; a = 2 describes to the mean-field case, whereas a = 1.33 is appropriate 
to the fluctuation dominated X — Y regime. This then yields the central prediction for the 
phase boundary B'^l'^{T) separating disordered liquid from multi-domain glass, [H] 

m(T — T 

Bt\T) = C{T - ^ - (5) 

In general, we could use any formula here which best models the melting line in the 
pure system, as inferred, for example, from fits close to Tc. Note that the suppression of 
the pure melting line by disorder is very weak if m is small. This suppression becomes 
progressively large as m is increased, or alternatively, at a lower temperature (larger H) 
for given m. It might appear that the final formula (Eq. [s]), with a, C and m as potential 
unknowns, contains a fairly large number of fitting parameters. However, note that as T — >■ 
Tc, the suppression term becomes irrelevant. Thus, two of the three parameters are fixed 
vis a vis the pure system, or equivalently, by an independent fit close to Tc; only a single 
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free parameter is required to fit the suppression of the melting hne by disorder, provided an 
independent fit to the pure system melting line close to B.c2 is available. 

We now describe the details of our experiments. The peak effect in pure 2H — NbSe2 
single crystals (Tc ~ 7. IK) is tracked by ac susceptibility (a 10 Hz ac field of amplitude 1 
Oe) and dc magnetization measurements using a Quantum Design SQUID magnetometer. 
Susceptibility measurements are carried out by varying temperature at different magnetic 
fields ranging from 50 Oe - 10000 Oe. Critical currents are inferred from field dependent 
magnetization hysteresis measurements at different temperatures above 2.5 K. The peak 
effect manifests itself in a susceptibility measurement through a dip in the real part of the 
ac susceptibility xu the minimum value of this quantity corresponds to the maximum value 
of the critical current. Typical plots of xi vs T at applied fields of 6000, 2000 and 200 
Oe is displayed in Fig. |3j Fig. |4] shows critical currents as inferred from magnetization 
measurements as the field is varied. The peak effect in magnetization measurements is 
obtained through the broadening of the magnetization hysteresis loops. The locus of the 
peak field vs temperature H(T) is independent of the technique used. In both cases a sharp 
peak is observed. 

Fig. |2] contains the main result of this paper, exhibiting points for the locus of peak 
positions {Hp,Tp) in H — T space, as obtained from susceptibility measurements. A fit 
to the freezing line for the multi-domain glass state, Eqn. [5} is shown in the same plot, 
with the fitting parameters C = 10^, m = 10.0, Tc = 7.09 and a = 1.33 (assuming XY 
exponents). The dashed line is a fit to the case where m = 0, as is conventionally assumed 
in Lindemann-parameter approaches to the melting transition. Observe the almost perfect 
agreement, over 3 decades or more in magnetic field JIBj. As we argue above, the suppression 
of the pure melting line can be parametrized by a single number, related to the suppression 
of the direct correlation function in the vicinity of the ordering wave-vector. Thus, our 
approach to parameterizing this phase boundary, apart from the agreement with experiment 
we demonstrate, is also economical in terms of being a single-parameter fit. 

In conclusion, we stress the main features of our analysis and our central result. We 
compute the suppression of the melting phase boundary in the pure system by quenched 
disorder. In contrast to approaches based on Lindemann parameter based measures of 
the instability of the Bragg glass phase, we determine the instability of the liquid to a 
static density wave. This instability, within the density functional formalism, arises as a 
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consequence of correlations which build up in the liquid phase; the method applied here 
parametrizes the effects of disorder in reducing such correlations. In general, it is to be 
expected that theoretical methods which study the instability towards freezing of a fluid 
phase should better equipped to capture the physics of the MG-DL line than methods 
which study the instability of a solid phase using a Lindemann criterion or its variants. 

Our central result, an expression for the disorder-suppressed phase boundary in the H — T 
plane across which the disordered liquid becomes unstable to a state which at least locally 
resembles a crystal, provides a remarkably accurate fit to the experimental data on the 
classic low-Tc superconductor 2H-NbSe2. The physical picture we outline is consistent with 
a general link between the loci of peak effect phenomena, reflecting the onset of complex 
dynamics of a driven vortex system, and underlying equilibrium order-disorder transitions as 
mirrored in the static phase diagram of Fig. [l |l3f[T5] . More work to test the validity of the 
theoretical prediction of B^^{T) against data from a variety of superconducting materials 
as well as to investigate the proposed universality of the phase diagram of Fig. [T] would be 
very useful. 

Support from DST (India) is gratefully acknowledged. We thank NEC Research Institute, 
Princeton, where these crystals were grown, for support. 
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